Strong saturation absorption imaging 
of dense clouds of ultracold atoms 
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We report on a far above saturation absorption imaging technique to investigate the characteristics of dense 
packets of ultracold atoms. The transparency of the cloud is controlled by the incident light intensity as a result 
of the non-linear response of the atoms to the probe beam. We detail our experimental procedure to calibrate 
the imaging system for reliable quantitative measurements, and demonstrate the use of this technique to 
extract the profile and its spatial extent of an optically thick atomic cloud. © 2008 Optical Society of America 

OCIS codes: 110.0110, 020.7010 
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Recently there has been a resurgent interest in the pro- 
duction of dense samples [1,2] containing a large number 
of cold neutral atoms with highly compressed magneto- 
optical traps (MOT) [3,4]. These studies, combined with 
optical trapping, opened the way to a simplified and very 
rapid production of Bose-Einstein condensates [5], and 
may also play a key role in the production of a cw atom 
laser based on the periodical coupling of atomic pack- 
ets into a magnetic guide, yielding a promising starting 
point for evaporative cooling [6]. 

For dense clouds, an important issue is the reliability 
of the method used to extract the atomic densities. The 
predominant imaging techniques for dilute samples, low- 
intensity fluorescent and absorption imaging, turn out to 
be unreliable when probing a dense atomic packet [7, 8] . 
The former critically depends on the value of the illu- 
minating intensity, the frequency and the repartition of 
atoms between the different Zeeman sub-levels. The lat- 
ter poses a problem as soon as the optical depth, propor- 
tional to the column density of atoms along the probe 
direction, is on the order of 3 to 4, because of both (i) 
the electronic noise of each pixel of the charge-coupled- 
device (CCD) camera and (ii) the digitalization of the 
signal of the weak amount of light that remains after the 
propagation through the cloud. 

An optically thick cloud can also be probed with off- 
resonant light or by implementing the phase contrast 
imaging technique [7]. The reduction of the scattering 
cross section that results from the non-resonant probe is 
favorable. However, the sample becomes dispersive and 
behaves like a gradient-index lens. Such a technique is 
fruitful for qualitative or differential measurements, but 
is difficult to use for quantitative purposes. 

To circumvent those drawbacks the group of D. Weiss 
has successfully implemented a fluorescence imaging 
technique very far above saturation intensity [3]: the 
probe intensity was larger by more than three orders of 
magnitude than the saturation intensity. In this regime, 



all atoms, independently of their Zeeman sub-level dis- 
tribution, spend half of the time in the excited state. 
This method allows to probe a cloud with a very large 
optical depth, but requires a dramatically high incident 
intensity and a subtle alignment of the two counter prop- 
agating beam used to drive the fluorescence. 

In this article, we report on our realization of a ro- 
bust, accurate and reliable far above saturation intensity 
absorption imaging aimed at investigating such dense 
atomic samples. Its implementation is straightforward 
and, by contrast with the far above saturation fluores- 
cence imaging technique, the probing does not required 
the use of a powerful laser. 

Indeed, in our experiment, the probe light is pro- 
vided only by diode lasers. A semiconductor slave laser 
is injection-locked to a 0.5 MHz linewidth Distributed 
Bragg Reflector (DBR) master diode laser, and spatially 
filtered by a pinhole. The purpose of this arrangement is 
to benefit from a narrow linewidth probe tuned on the 
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tively large power (30 mW) available to probe the atoms. 
An acousto-optic modulator placed before the pinhole 
is used to produce light pulses as short as 250 ns. The 
shadow cast by the atoms on the resonant probe beam 
is imaged on a CCD camera, with an optical resolution 
of 7 i-j.ni. 

The response of the atoms, i.e. the population driven 
in the excited state by the imaging laser beam, depends 
on the effective saturation intensity J|g = a*/o a \ where 
/q 11 ' is the saturation intensity for the corresponding two- 
level transition (Ig at = 1.67 mW/cm 2 for rubidium 87). 
The dimensionless parameter a* accounts for corrections 
due to the specific conditions in which images are taken: 
the polarization of the imaging beam, the structure of 
the excited state and the different Zeeman sub-level pop- 
ulations of the degenerate ground state of the optical 
transition. 

In order to extract the spatial atomic density n(x, y, z) 
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of the cloud, wc acquire as usual three images: I w (x, y) 
with the atoms and probe beam on, I wo (x, y) without the 
atoms and probe on, and /dark (x, y) without atoms and 
probe off. From those images, we work out for each pixel 
(x, y), the light intensity if (a:, y) = I w (x, y) - /dark (x, y) 
(resp. h{x,y) = I wo (x,y) - Id a ik(x,y)) of the imaging 
beam in presence (resp. absence) of atoms by remov- 
ing the contribution of the background light illumination 
taken in the absence of the detection beam. 

The Beer's law in presence of saturation effect and for 
a resonant incident light can be recast in the form: 

dI g ° 1 t - it\t m 
— = — n . , 1 = —71(7(1)1 (1) 

dz a* l + I/I^ y ' w 

where <7q = 3\ 2 /2ir is the resonant cross-section for a 
two- level atom, and a (I) the effective cross section in- 
cluding saturation correction. From Eq. (1), one readily 
obtains the expression for the optical depth: 

od (x, y)=a J n(x, y, z) dz = f(x, y; a*). (2) 

where f(x,y;a*) is defined by: 

(3) 

The optical density is defined by: So(x,y) = 
— ln(If(x,y)/Ii(x,y)). We have deliberately chosen a 
definition of the optical depth odo(x,y) that does not 
include the parameter a* so that it depends only on the 
characteristics of the atomic cloud. However, to extract 
this quantity from a set of images, one needs to know 
the parameter a*. 

For low- intensity absorption imaging (I[(x, y) <C io at ), 
the optical depth involves only the ratio of the intensi- 
ties If and /: odo(x,y) ~ a*6o(x, y). The unknown pa- 
rameter a* still needs to be determined independently. 
The strong saturation imaging technique takes advan- 
tage of the reduction of the effective cross-section a (I) 
when I 3> a*lQ at . In this limit, the optical depth de- 
pends both on a* and on the value of the incident inten- 
sity Ii(x,y), because of the non- linear atomic response 
(see Eq. (3)). 

The absolute calibration is a crucial but delicate task 
for fluorescence or low absorption imaging. In the con- 
text of our far above saturation absorption imaging tech- 
nique, it just consists in determining the parameter a*. 
We proceed in the following manner. The sample of 
cold atoms is generated by a compressed elongated two- 
dimensional MOT. The cloud is imaged after a not too 
short time-of-flight so that its maximum optical den- 
sity is not too high (~ 2), which guarantees also the 
validity of the low-intensity absorption imaging. As the 
imaging technique is destructive, we acquire several set 
of images (typically five) of a cloud always prepared in 
the same conditions for different incident intensities. In 
practice, we vary the intensity of the imaging beam by 
more than two orders of magnitude while keeping the 



number of photons per pulse constant: the duration of 
the pulse was varied from 250 ns (/ ~ 23 mW/cm 2 ) to 
100 fis (Ii ~ 0.06 mW/cm 2 ). Keeping the number of ab- 
sorbed photons small (~ 5 photons per atom on average) 
avoids pushing and heating the cloud, thereby changing 
its characteristics. 

In order to infer the value of the dimensionless param- 
eter a* , we calculate, for different values of the parame- 
ter a ranging from 1 to 4, the function f(x, y; a) for the 
set of images. We extract the amplitude od(a) of those 
calculated optical depth using a gaussian fit (see Fig. 1). 
There is only one value a* of a for which all the calcu- 
lated od(a) are equal over the whole range of incident 
intensities. Indeed, odo does not depend on the incident 
probe intensity (see Eq. 2). 

In practice, we infer by a least square method the value 
a* for which od(a) has a minimum standard deviation 
over the whole range (more than two orders of magni- 
tude) of incident intensities used to image the cloud (see 
Fig. 1 inset). We find a* = 2.12 ± 0.1 and a maximum 
optical depth od = 4.8 that corresponds to an opti- 
cal density of <5 = 2.25 as deduced from the low inten- 
sity absorption imaging. We stress that this calibration 
allows for an absolute determination of the number of 
atoms, its accuracy being ultimately determined by the 
calibration of the incident intensity using the CCD array. 

This value for a* is to be compared with the result 
of the Bloch equations for the corresponding multiple 
level system. For our data, the atoms arc initially in 
the \g) = 5 2 Si/2,F = 2 hyperfine state (5-fold degen- 
erated) and the 7r-polarized probe is resonant with the 
|e) = 5 2 P 3 / 2 ,F' = 3 hyperfine excited state (7-fold de- 
generated). The probability to excite the 5 2 -P 3 / 2 , F' = 2 
is negligible (below 1 %). The transition \g) — > |e) can 
be considered in this limit as closed. From numerical in- 
tegration of the Bloch equations, we find out that (i) 
the steady state solution is approximately valid even for 
the shortest pulses that we use, and (ii) the correction 
factor a* to the two-level saturation intensity lies in be- 
tween one and two depending on the polarization of the 
probe. Experimentally, the polarization of the beam can- 
not be perfectly under control because of the slight bire- 
fringence of the viewports. In addition, a residual mag- 
netic field may also influence the effective value of a* . 

After calibration, the high-intensity imaging technique 
is applied to the same cloud but without time-of-flight. 
The transparency of the atomic cloud is controlled by the 
probe intensity because of the non-linear response of the 
atoms. Our numerical studies based on the Eqs. (2) and 
(3) show that the true profile of the cloud can be reliably 
inferred from the strong saturation absorption images as 
soon as the incident intensity of the probe beam is on the 
order of / ~ max(o<io) 1^. Note that this intensity is 
much less than the one needed for a reliable high inten- 
sity fluorescence imaging technique [3] . We have checked 
this prediction by imaging with different incident inten- 
sities a compressed two-dimensional MOT prepared in 
the same conditions. To compress the MOT, we proceed 
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in the following manner: the repumper intensity is di- 
vided by 15 in 1 ms, the detuning is ramped linearly 
from — 3r to — 9r in 15 ms, and the gradient is increased 
from 5 G/cm to 20 G/cm in 15 ms. From our analysis 
based on the strong saturation absorption, we system- 
atically find out a double structure (see Fig. 2.c) with 
a central dense region having an maximum peak atomic 
density on the order of 2(±1) x 10 11 atoms/cm 3 corre- 
sponding to max(odo) ~ 45. Half of the atoms (1.5 x 10 8 ) 
remains in the wing. The profile of those wings and their 
number of atoms inferred from the low and high intensity 
absorption imaging techniques perfectly coincides. 

In addition, this technique is particularly well suited 
for the estimation of the spatial extent of dense atomic 
packets with one or more sizes very small (Axq < 30 
//m). If T is the temperature of the cloud, the size of the 
packet reflects the velocity distribution after a time-of- 
flight duration on the order of t = (mAa^/ZcsT) 1 / 2 . For 
a cloud of initial size 10 /Ltm and a temperature of 100 
/xK, r ~ 100 fj,s. By contrast with low-intensity absorp- 
tion imaging, for which the pulse duration is on the order 
of the time t, the very short pulse used for high-intensity 
imaging prevents heating and permits to extract the cor- 
rect profile. 

We have demonstrated a high-intensity absorption 
imaging technique well-suited for dense, small-size 
atomic clouds. This technique is robust against small fre- 
quency and intensity variations of the probe, and does 
not require an extremely large probe intensity. We have 
shown how it can be reliably calibrated exploiting the 
non-linear atomic response. 
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Fig. 1. (color online) A cloud is imaged using differ- 
ent probe intensities (from /Q at /15 to 15/o at ). For each 
image, the maximum optical depth of the cloud od(a), 
deduced from a gaussian fit, is calculated with several 
values of the unknown parameter a using the function 
/. The plot represents od(a) as a function of the incom- 
ing intensity. The standard deviation A(od) of each set 
of data points (see inset) exhibits a clear minimum as 
a function of the parameter a. The minimum of A(od) 
gives a* = 2.12 ±0.1. 
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Fig. 2. Three images with a transverse cut of a dense 
elongated cloud prepared in the same conditions. The 
transverse profile is given by the dimension less func- 
tion f(x,y;a*)/a* and the dashed line corresponds to 
the value 3. (a) Low intensity absorption imaging. Al- 
most all the light is absorbed in the center of the cloud, 
(b) Off-resonance absorption imaging with a probe laser 
blue-detuned by 1.2r where T is the natural linewidth of 
the excited state. One observes a double peak structure 
resulting from a gradient-index lens effect, (c) High in- 
tensity absorption imaging. Only the third method shows 
that the profile has a double structure, with an r.m.s 
size for the central peak of 300 /xm, corresponding to 
a peak atomic density ~ 4(±1) x 10 10 atoms/cm 3 and 
max(odo) ~ 9. 
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